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1.  Introduction 


US  Army  Research  Laboratory  (ARL)  researchers  have  achieved  robust  theoretical 
models  capable  of  predicting  performance  properties,  such  as  heats  of  formation1,2 
and  densities3,4  of  energetic  materials.  This  capability,  coupled  with  a  growing  ARL 
synthetic  experimental  faculty,  allows  synthetic  and  formulation  chemists  to  safely 
and  quickly  screen  candidate  materials  to  focus  efforts  only  on  the  most  promising 
compounds.  This  technical  note  will  detail  theoretical  predictions  of  a  notional 
energetic  material. 

The  tools  are  based  on  correlations  of  quantum  mechanically  derived  properties  of 
isolated  energetic  molecules  with  heats  of  formation  and  crystalline  densities, 
critical  elements  in  predicting  the  detonation  and  propulsion  properties,  for  both 
neutral  and  ionic  solid  energetic  materials.  As  the  crystal  structure  of  a  novel 
material  is  usually  unknown,  as  it  has  yet  to  be  synthesized,  each  of  the  correlations 
requires  quantum  mechanical  data  solely  from  single  molecules  and  not  from 
knowledge  of  the  bulk  material.  To  date,  the  neutral  heat  of  formation  tool  can  only 
treat  molecules  containing  carbon,  hydrogen,  nitrogen,  and  oxygen.  Additionally, 
a  tool  to  qualitatively  predict  the  impact  sensitivity  for  neutral  molecules  was 
developed,  also  based  solely  on  quantum  mechanically  derived  data  based  on  a 
single  molecule.  This  report  details  the  theoretically  predicted  performance  and 
vulnerability  properties  for  notional  energetic  molecules;  as  such,  we  will  recap  the 
methodologies  used  in  the  predictions  in  this  technical  report. 

1.1  Neutral  Heat  of  Formation 

Developed  in  2006, 1  the  neutral  heat  of  formation  prediction  tool  is  based  on  Hess’s 
law  (Eq.  I),5  where  an  atomic  or  group  equivalency  prediction  for  the  gas  phase 
heat  of  formation  (Eq.  2)  is  added  to  a  prediction  for  the  heats  of  vaporization 
(Eq.  3)  or  sublimation  (Eq.  4).  Equations  2-4  use  information  determined  from 
quantum  mechanical  calculations. 


AH°f  (liquid/sold)  —  AH°f(g)  -  AHvap/sub- 

(1) 

AH°f(g)  =  E  -  I  rijEj. 

(2) 

AHvap  =  a(SA)1/2  +  b(a2to,v ) 1/2  +  c. 

(3) 

AHsub  =  a(SA)2  +  b(o2totV ) 1/2  +  c. 

(4) 

In  Eq.  2,  E  is  the  quantum  mechanically  determined  electronic  energy,  nj  is  either 
the  number  of  atom  or  group  types  j  contained  in  the  neutral  molecule,  and  sj  is 
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either  the  atom  or  group  equivalent  energy  of  atom/group  j.  These  equivalent 
energies  were  determined  through  the  fitting  of  experiment  gas  phase  heats  of 
formation  [AH°f  (g)]  to  quantum  mechanically  computed  energies  E. 

Equations  3  and  4  use  statistics  derived  from  the  electrostatic  surface  potential 
(ESP),  as  demonstrated  by  Politzer  et  al.6-8,  to  correlate  to  heats  of  vaporization  and 
sublimation.  As  defined  in  the  Murray  et  al.9  and  Rice  and  Hare10  papers,  the 
definitions  of  these  electrostatic  properties  are  global  properties  across  the  0.001 
electron/bohr3  isosurface  of  electronic  density,  an  isosurface  calculated  to  contain 
98%-99%  of  the  electronic  density.  To  compute  these  statistics,  the  electrostatic 
potential  V  is  required  and  is  defined  as 


v(r)=Y 

v  }  L—1[  lR;-rl  J  Ir'-rl  ’ 


(5) 


where  Z,  and  Ri  represent  the  charge  and  position  of  atom  i,  and  p{ r)  is  the  electronic 
density.  The  statistical  tools  developed  by  Murray  et  al.9,11  used  in  Eqs.  3  and  4  are 
the  surface  area  (SA),  the  variance  of  electrostatic  surface  potential  (a2),  and  the 
balance  parameter  (v).  These  are  defined  as 

•  (6) 
°l  =  ,  (7) 

=  a't+  ai  ,  (8) 


and 


v  - 


2  ? 


(9) 


where  V(n)  is  the  value  of  the  ESP  at  r,  and  X  is  the  average  electrostatic  surface 
potential  (x  is  either  the  positive  [+]  or  negative  [-]  portion  of  the  ESP).  SA  is  the 
surface  area  of  the  0.001  electron/bohr3  isosurface  of  the  electron  density  of  the 
molecule,  o2t0t  is  a  measure  of  the  variability  of  electronic  potential  on  the  surface, 
and  v  is  the  degree  of  balance  between  the  positive  and  negative  charges  on  the 
isosurface.  Equations  3  and  4  contain  a  different  set  of  fitted  a,  b,  and  c  parameters. 
This  methodology  enjoys  considerable  accuracy,  with  the  root  mean  square  errors 
for  the  solid  phase  heat  of  fonnation  under  6  kcal/mol. 
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1.2  Crystalline  Density 


Initially  developed  in  20073  and  revised  in  2013, 4  this  method  correlates  the 
quantum  mechanically  derived  electronic  density  volume  of  the  molecular  unit  to 
the  crystalline  density.  This  correlation  exhibits  significant  success  in  predicting 
the  crystal  density,  with  errors  under  3%  for  neutral  materials.  This  correlation  only 
requires  the  quantum  mechanically  derived  molecular  volume  and  the  molecular 
mass  for  a  single  molecular  unit,  along  with  associated  electrostatic  properties 
(Eq.  10). 


P  crystal  a(M  /  Vm)  +  b(vo2tot )  +  c,  (10) 

where  a,  b,  and  c  are  fitted  terms,  different  from  those  in  Eqs.  3  or  4,  M  is  the  total 
molecular  mass,  Vm  is  the  total  molecular  volume,  and  o2t0t  and  v  are  the  same  as 
in  Eqs.  8  and  9. 

1.3  Impact  Sensitivity 

Besides  perfonnance  parameters,  ARL  researchers  have  also  developed  a  visual 
qualitative  indicator  of  impact  sensitivity  (David  Chavez,  Los  Alamos  National 
Laboratory,  private  communication,  2014;  unreferenced)  as  a  method  to  assess  any 
potential  safety  concerns  with  yet  to  be  synthesized  novel  materials.  Rice  and 
Hare12  attempted  to  correlate  Politzer-type  statistical  properties  of  the  electronic 
density  to  impact  sensitivity  for  a  series  of  polynitroaromatic  and  benzofuroxan 
molecules.  While  no  strong  correlations  were  found  with  the  statistical  parameters, 
there  was  a  qualitative  visual  correlation  related  to  the  degree  of  electron  deficiency 
over  covalent  bonds  contained  within  the  “backbone”  of  the  molecules.  Molecules 
that  were  very  sensitive  exhibited  large  electron  deficiencies  over  the  main 
structure,  while  less-sensitive  molecules  had  a  more  evenly  distributed  electron 
density.  This  method,  while  quite  approximate  and  limited  to  neutral  molecules 
(due  to  the  inherent  charge  separation  found  in  ionic  materials),  can  still  yield  a 
relative  means  of  estimating  impact  sensitivities  and  thus  potential  hazards  of 
notional  materials.  The  visualizer  used  in  this  study,  titled  ED  AT,  was  originally 
coded  by  Dr  Jennifer  Hare  and  later  revised  by  Dr  Anthony  Yau.  This  tool  is  also 
used  to  derive  the  Politzer  statistical  properties  used  in  the  heat  of  formation  and 
density  predictors.  In  this  report  we  shall  apply  these  methods  toward  the  molecule 
3,7-dinitro-[l,2,4]triazolo[5,l-c][l,2,4]triazin-4-amine  (David  Chavez,  Los 
Alamos  National  Laboratory,  private  communication,  2014;  unreferenced)  to 
predict  the  heat  of  formation  and  crystalline  density  of  this  material. 
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2.  Method  of  Solution 


The  analysis  was  performed  using  scripts  written  at  ARL  to  perform  the  following 
required  computations12: 

1 .  A  driver  script  that  calls  the  other  scripts,  as  well  as  submits  the  required 
Gaussian0913  calculations  and  monitors  the  queues  (master  densheat.s.) 

2.  A  script  to  construct  the  Gaussian  input  fdes  (gbuild.s). 

3.  A  script  to  gather  and  compute  the  data  needed  to  determine  the  heats  of 
formation  and/or  crystalline  densities  and  subsequently  calculate  said 
properties  (arl  densheat.s). 

Given  a  properly  fonnatted  input  file,  the  scripts  will  analyze  the  .dat  input  file, 
design  the  proper  Gaussian  input  files  and  submit  them,  monitor  the  supercomputer 
queues  for  job  completion,  postprocess  the  necessary  files,  and,  finally,  compute 
the  desired  information. 

The  scripts  also  perform  a  limited  series  of  error  checking,  such  as  checking  that 
structures  are  optimized  (i.e.,  the  structure  is  not  a  transition  state  structure  but 
rather  a  minimum),  that  the  Gaussian  calculations  ran  successfully  to  completion, 
and  that,  for  neutral  calculations,  the  number  of  atom  types  listed  correspond  to  the 
number  of  atoms  in  the  Cartesian  coordinates. 

The  initial,  user-designed  structure  has  its  geometry  optimized  using  the  B3LYP 
spin-restricted  Kohn-Sham  density  functional  theory  (KS-DFT)14  with  the  6-3 1 G* 
Pople  Gaussian  basis  set,15  currently  done  using  the  Gaussian09  program 
package.13  For  the  computation  of  the  electronic  energy  required  in  Eq.  2,  the 
B3LYP  density  functional  theory  energy15  with  the  6-31 1++G(2df,2p)  Pople  basis 
set16  is  obtained  on  the  optimized  geometry.  For  the  crystalline  density,  the  initial 
structure  is  optimized  using  the  B3LYP  KS-DFT  with  the  6-3 1G**  Pople  basis 
set,16  with  the  associated  electronic  density  used  in  obtaining  the  molecular  volume 
and  electronic  statistical  properties  required  in  Eq.  10.  The  visualization  of  the 
electrostatic  potential  surface  was  performed  with  the  ED  AT  tool. 

3.  Results  and  Discussion 


The  3,7-dinitro-[l,2,4]triazolo[5,l-c][l,2,4]triazin-4-amine  molecule  (David 
Chavez,  Los  Alamos  National  Laboratory,  private  communication,  2014; 
unreferenced)  (Fig.  1)  was  constucted  and  submited  with  the  ARL  scripts.12 
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Fig.  1  Optimized  structure  of  the  3,7-dinitro-[l,2,4]triazolo[5,l-c][l,2,4]triazin-4-amine 
molecule 

The  resultant  crystal  density  was  computed  to  be  1.856  g/cc.  The  gas  and  solid 
phase  heats  of  formation  from  the  group  additivity  methodology  were  120.12  and 
95.22  kcal/mol,  respectively. 

Using  the  ARL-developed  ED  AT  tool,  the  electrostatic  map  on  the  0.001  isosurface 
was  generated  and  is  illustrated  in  Figs.  2  and  3  (with  and  without  the  molecule 
overlaid).  For  this  visualization,  the  scalar  range  of  the  ESP  ranges  from  -0.05  to 
0.075.  For  this  visualization  methodology,  regions  of  large  positive  charge  (red) 
over  the  backbone  of  the  structure  tend  to  indicate  increased  sensitivity.  For  this 
molecule,  this  ESP  distribution  would  indicate  moderate  impact  sensitivity.  Using 
this  predicted  data,  we  then  performed  Cheetah  7.0  calculations18  to  predict  the 
performance  parameters  (Table). 
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Fig.  2  Electrostatic  potential  map  of  the  3,7-dinitro-[l,2,4]triazolo|5,l-c][l,2,4|triazin-4- 
amine  molecule,  with  molecule  overlay 


Fig.  3  Electrostatic  potential  map  of  the  3,7-dinitro-[l,2,4]triazolo|5,l-c][l,2,4|triazin-4- 
amine  molecule,  without  molecule  overlay 
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Table  Cheetah-predicted  properties  of  the  3,7-dinitro-[l,2,4]triazolo|5,l-c]|l,2,4]triazin-4- 
amine  molecule 


Property 

Value 

Pressure 

31.084  GPa 

Shock  velocity 

8.679  km/s 

Temperature 

Total  energy  of  detonation 

3,395.8  K 

1.284  TNT  equivalence  (per  cc) 

1.145  TNT  equivalence  (per  gram) 

4.  Conclusions 


The  ARL-developed  software  tools  were  used  to  predict  the  heat  of  formation  and 
crystalline  density  of  the  3,7-dinitro-[l,2,4]triazolo[5,l-c][l,2,4]triazin-4-amine 
molecule.  Using  this  predicted  data,  we  then  ran  Cheetah  calculations  to  predict  the 
performance  of  this  novel  material.  Additionally,  we  predicted  the  qualitative 
impact  sensitivity  of  this  compound  using  electrostatic  potential  maps.  This 
information  has  been  transferred  back  to  the  requesting  synthetic  chemist,  Dr  David 
Chavez  of  Los  Alamos  National  Laboratory. 
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List  of  Symbols,  Abbreviations,  and  Acronyms 


6-3 1G* 

Pople’s  Gaussian  double-zeta  polarized  basis  set  with  d 
polarization  functions  on  each  of  the  atoms  Li  through  Ca 

6-3 1G** 

Pople’s  Gaussian  double-zeta  polarized  basis  set  with  d 
polarization  functions  on  each  of  the  atoms  Li  through  Ca 
and  p  polarization  functions  on  H  and  Fie 

6-31 1++G(2df,2p) 

Pople’s  Gaussian  triple-zeta  polarized  basis  set  with  2d  and 
/ polarization  functions  on  each  of  the  atoms  Li  through  Ca 
and  2 p  polarization  functions  on  H  and  He  along  with  s  and 
p  diffuse  functions  on  each  of  the  atoms  Li  through  F  and  s 
diffuse  functions  on  H  and  He 

V 

balance  parameter 

I 

average  electrostatic  surface  potential 

a2 

variance  of  electrostatic  surface  potential 

P 

electronic  density 

Pcrystal 

crystalline  density 

AH°f  (g) 

gas  phase  heat  of  formation 

AH°f(l) 

liquid  phase  heat  of  formation 

AH°f(s) 

solid  phase  heat  of  fonnation 

V 

electrostatic  potential 

ARL 

US  Army  Research  Laboratory 

B3LYP 

Becke  3  parameter  exchange  with  Lee- Yang-Parr 
correlation  DFT  functional 

DFT 

density  functional  theory 

ESP 

electrostatic  surface  potential 

g/cc 

grams  per  cubic  centimeter 

kcal/mol 

kilocalories  per  mole  (unit  of  energy) 

KS-DFT 

Kohn-Sham  density  functional  theory 
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